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Abstract 

A  list  of  publications  of  the  research  performed  during  the  period  1/1/80- 
6/30/84  is  provided.  Theoretical  research  has  been  conducted  (a)  on  ion-ion 
recombination  in  high  pressure  plasmas,  (b)  on  the  theory  of  Rydberg 
collisions  with  electrons,  ions  and  neutrals,  and  (c)  on  the  theori  of  ion- 
molecule  collisions  at  (1  eV-4  keV)/AHU.  Papers  in  these  topics  (a)-(c)  have 
been  written  up  and  published  as  papers  already  submitted  to  AFOSR.  In 
particular,  a  new  and  basic  theory  of  ion-ion  recombination  in  a  dense  gas  has 
been  developed  from  basic  microscopic  principles.  The  recombination  rate  is 
provided  as  a  function  of  gas  density  and  time.  The  theory  has  been  published 
in  papers  already  submitted  to  AFOSR.  ^pendix  A  of  this  report  contains  the 
published  chapter  on  "Theory  of  Ion-Molecule  Collisions  at  (1  eV-50  keV)/AMU". 
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M.  R.  FLANNERY 

School  of  Physics,  Georgia  Institute  of  Technology 

Grant  AFOSR-80-0055 

Project  No.  2301,  Task  No.  A4 


1.1  Full  List  of  Refereed  Publications  In  Scientific  Journals  (1980-1984) 

1.  "Charge-Transfer  In  Three-Body  Ion-Ion  Recombination  at  Low  Gas 
Densities",  Int.  J.  Quant.  Chem.:  Quant.  Chem.  Symp.  14,  477-482  (1980). 

2.  "Ion-Ion  Recombination  In  Systems  at  Low  Gas  Densities:  I. 

Symmetrical  Resonance  Charge-Transfer  Contribution",  J.  Phys.  B:  Atom. 
Molec.  Phys.  13,  3649-3664  (1980). 

3.  “Three-Body  Recombination  of  Rare-Gas  Atomic  Ions  X'*’  with  F"  In  a 
Low-Density  Gas  X",  with  T.  P.  Yang,  J.  Chem.  Phys.  73,  3239-3245  (1980). 

4.  "Theoretical  Treatment  of  Collisions  of  Rydberg  Atoms  with  Neutral  Atoms 
and  Molecules.  The  Semiquantal,  Impulse  and  Multi state-Orbltal 
Theories",  Phys.  Rev.  A  2408-2429  (1980). 

5.  "Thermal  Collisions  of  Rydberg  Atoms  with  Neutrals",  J.  Phys.  B:  Atom. 
Molec.  Phys.,  L657-663  (1980). 

6.  "Ion-Ion  Recombination  In  (X^+Y'+Z)  Systems  at  Low  Gas  Densities.  II. 
Elastic 

Ion-Neutral  Collisions",  J.  Phys.  B:  Atom.  Molec.  Phys.  915-934 
(1981). 

7.  "Vibrational  Deactivation  of  Oxygen  Ions  In  Low  Velocity  02^(X  v“l) 

+  02(X'*  e",  v“0)  Collisions",  with  T.  F.  Moran,  K.  J.  McCann,  M.  Cobb  and 
R.  F.  BorRman,  J.  Chem.  Phys.  74,  2325-2330  (1981). 

8.  "Ion-Ion  Recombination  as  a  Function  of  Ion  and  Gas  Densities",  Chem. 
Phys.  Letts.  M,  541-546  (1981). 

9.  "Exact  Closed  Form  Solution  of  the  Generalized  Debye-Smo1uchowsk1 
Equation",  Phys.  Rev.  Letts.  47,  163-166  (1981). 

10.  "Cross  Sections  for^Elect^on-Impact  Ionization  of  Metastable  Rare-Gas 
Excimers  (He,  ,  Kr,  ,  Xe-  )",  with  K.  J.  McCann  and  N.  W.  Winter,  J. 

Phys.  B:  Atom.  Mofec.  Phys.  H,  3789-3796  (1981). 

11.  "Ion-Ion  Recombination  in  Dilute  and  Dense  Plasmas",  Int.  J.  Quant. 

Chem.:  Quant.  Chem.  Symp.:  715-727  (1981). 

12.  "Theory  of  Ion-Ion  Recombination",  Phil.  Trans.  Royal  Soc.  (London)  A 


304,  447-497  (1982) 


13.  "Ion-Ion  Recombination  In  High  Pressure  Plasmas",  In  Applied  Atomic 
Collision  Physics.  Vo1«  3:  Sas  Lasers.  Eds.  H.  S.  W.  Massey,  B. 

Bederson  and  E.  W.  McDaniel,  (New  Vork,  Academic,  1982)  Chapter  5,  pages 
141-172. 

14.  "Analytical  Solutions  of  the  Debye-SmoluchoMSkl  Equation  for  Geminate  and 
Homogeneous  Recombination  and  for  Fluorescence  Quenching",  Phys.  Rev.  A. 
K,  3403-3406  (1982). 

15.  "Comments  on  Thermal  Collisions  of  Rydberg  Atoms  with  Neutrals",  J.  Phys. 
B:  Atom.  Molec.  Phys.  3249-3256  (1982). 

16.  "Exact  Closed  Form  Solution  of  the  Generalized  Debye-Smoluchowski 
Equation",  ^ys.  Rev.  Letts.  49,  1681  (1982). 

17.  "Theory  of  Rydberg  Collisions  with  Electrons,  Ions  and  Neutrals",  In 
Rydberg  States,  Eds.:  R.  F.  Stebbings  and  F.  B.  Dunning  (Cambridge, 
Cambridge  University  Press,  1983),  Chapter  11,  pp.  393-454. 

18.  "Diffuslonal-Orlft  Influenced  Chemical  Reactions  at  all  Gas  Densities", 
Phys.  Rev.  Letts.  1656-1659  (1983). 

19.  "Theory  of  Ion-Molecule  Collisions  at  (1  eV-5  keV)/AMU",  In  Chapter  6  of 
Swarms  oF  Ions  and  Elytrons  In  Gases,  1984,  eds.  VI.  Lindinger,  T.  D. 
Mark  and  F.  Howorka  (Springer-Verlag),  pages  103-125. 

During  the  period  (1980-1984)  of  the  Grant  AFOSR-80-0055,  six  copies  of 
reports  describing  the  above  published  research,  II  -  118,  were  submitted  to 
AFOSR  routinely  as  they  became  available.  Also  full  details  were  contained  In 
the  Interim  Annual  Technical  Reports  submitted  to  AFOSR  for  the  periods  (a) 
January  1,  1980  -  Oecember  31,  1980,  (b)  January  1,  1981  -  December  31,  1981 
and  (c)  January  1,  1982  -  Oecember  31,  1982. 

The  published  article  #19  above  entitled  "Theory  of  Ion-Molecule 
Collisions  at  (1  eV  -  5  keV)/AMU",  Is  Included  as  Appendix  A  of  this  Final 
Report. 
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.  "Ion-Ion  Recombination  in  High  Pressure  Plasmas",  Chapter  5  of  Applied 
Atomic  Collision  Physics.  Vol.  3,  Gas  lasers.  1982  eds.,  H.  S.  W.  Hassey, 
B.  Bederson  and  E.  W.  McDaniel  (New  York,  Academic),  pages  141-172. 

2.  "Theory  of  Rydberg  Collisions  with  Electrons,  Ions  and  Neutrals",  in 
Chapter  II  of  Excited  Rydberg  States.  1983  eds,,  R.  F,  Stebbings  and  F. 

B.  Dunning  (Cambridge  University  Press),  pages  393-453. 

3.  "Theory  of  Ion-Molecule  Collisions  at  (1  eV-5  lteV)/AMU",  in  Chapter  6  of 
Swarms  of  Ions  and  Electrons  in  Gases.  1984,  eds.  W.  Lindinger,  T.  0. 
MarK  and  F.  Howorka  (Sprlnger-Veriag),  pages  103-125. 

These  chapters  have  also  been  included  in  the  list  (1.1)  of  publications 
but  are  identified  explicitly  here.  Six  copies  of  reprints  of  above  articles 
#1  and  #2  have  been  submitted  to  AFOSR  when  they  became  available.  Other 
details  were  contained  in  the  Interim  Annual  Reports  for  the  periods  (a) 
January  1,  1980  -  Oecember  31,  1980,  (b)  January  1,  1981  -  December  31,  1981 
and  (c)  January  1,  1982  -  Oecember  31,  1982.  Chapter  #3  is  included  as 
Appendix  A  of  this  report. 


1.3 .  Funding  History 

(a)  Total  Number  of  Publications  due  to  AFOSR  support  (1980-1984):  19 


(b)  Funding  FY  80  ; 

FY  81  ; 

FY  82  ; 

FY  84 

(1/1/80-12/31/80); 

(1/1/81-12/31/81); 

(1/1/82-12/31/82); 

(1/1/83-6/30/84) 

$51,000  ; 

$58,266 

$64,024  ; 

$69,919 

2.  Invited  Papers  Presented  at  Conferences  (1974-1984) 

In  addition  to  presentation  of  many  contributed  papers  at  the  various 

Annual  Meetings  (sponsored  by  the  American  Physical  Society)  of  the  Gaseous 

Electronics  Conference  and  The  Division  of  Electron  and  Atomic  Physics,  the 

following  Invited  Papers  were  presented  on  research  due  to  AFOSR  support. 

1.  “Atomic  and  Molecular  Collision  Processes  In  High-Power  Excimer  Lasers", 
Invited  Lecture  presented  at  International  Symposium  on  Atomic  Molecular 
and  Solid  State  Theory.  Florida,  March  11-  17.  1979. 

2.  "Theory  of  Excimer  Lasers",  Invited  Lecture  for  New  Directions  In  Atomici 
Physics,  Gordon  Research  Conference,  Brewster  Acaden^,  New  Hampshire,  June 
25-29,  1980. 

3.  "Basic  Microscopic  Theory  of  Ion-Ion  Recombination",  Invited  Lecture 
presented  at  International  Symposium  on  Atomic,  Molecular  and  Solid  State 
Theory,  Florida,  March  9-14,  1981. 

4.  "Basic  Microscopic  Theory  of  Neutralization  and  of  Chemical  Reactions  In 
Dense  Gases",  Special  Long  Paper  (30  mins.)  delivered  at  34  Annual  Gaseous 
Electronics  Conference,  boston.  Mass.,  October  20-23,  1981. 

5.  "Basic  Microscopic  Theory  of  Chemical  Reactions  In  Dense  Gases",  Invited 
Lecture  delivered  at  the  3rd  Annual  Symposium  on  Atomic  and  Surface 
Plasmas,  Hlntermoos,  Austria,  February  7-13,  1982. 

6.  "Ion-Ion  Recombination  In  Dense  Gases",  Invited  Paper  delivered  to 
International  Conference  on  Lasers  *82,  New  Orleans,  La.,  December  13-17, 
198?: 

7.  "Ion-Ion  Recombination",  Invited  Lecture  in  the  Distinguished  Invited 
Lectureship  Series  of  Harvard  University,  Feb.  1983. 

8.  "Theory  of  Ion-Ion  Recombination",  Invited  Plenary  Lecture  delivered  to 
International  Symposium  on  Atomic,  Molecular  and  Solid  State  Theory, 
Florida,  March  3-6,  1983. 


.  "Theory  of  Ion-Molecule  Collisions",  Invited  Lecture  at  Workshop.  XVth 
Annual  Meeting  of  the  Division  of  Electron  and  Atomic  Physics  (American 
Physical  Society),  May  30  -  June  2,  1984. 
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3.  Special  Highlights;  New  Theoretical  Developments  In  Present  Research 
A  new  and  basic  theory  of  ion-ion  recombination  as  a  function  of  gas 
density  N  has  been  developed  (M.  R.  Flannery,  Phil.  Trans.  Roy.  Soc.  A  304, 
447-497  (1982))  from  basic  microscopic  principles.  A  key  equation  for  the 
distribution  in  phase  space  of  ion  pairs  had  been  derived  together  with  an 
expression  for  the  resulting  recombination  coefficient  a.  Further  development 
of  the  theory  led  to  interesting  insights  to  a  full  variation  with  N  of  a, 
which  was  shown  to  yield  the  correct  limits  at  low  and  high  N.  The 
recombination  rate  a  is  determined  by  the  limiting  step  of  the  rate  for  ion 
reaction  and  of  the  rate  for  ion  transport  to  the  reaction  zone.  An 
accurate  analytical  solution  of  the  time-dependent  Debye-Smolcuhowski  equation 
which  is  a  natural  consequence  of  this  theory,  was  provided,  for  the  first 
time,  for  transport/reaction  under  a  general  interaction  V  in  the  cases  of  an 
instantaneous  reaction  »  o^^^)  and  of  a  finite  rate  reaction 

within  a  kinetic  sink  rendered  compressible  via  variation  of  gas  density. 
Expressions  for  the  transient  recombination  rates  ®(t)  were  then  derived  and 
illustrated.  The  exhibited  time  dependence  lends  itself  to  eventual 
experimental  verification  at  high  N. 

A  theory  which  investigated  the  variation  of  o  with  ion  density  was 
also  developed.  Here  the  ion-ion  interacton  V  can  no  longer  be  assumed 
ab-initio  to  be  pure  Coulomb  but  is  solved  self-consistently  with  the 
recombination.  Recombination  rates  for  various  systems  were  illustrated  as  a 
function  of  N  via  a  simplified  method  for  the  reaction  rate.  Finally,  two 
theoretical  procedures  for  the  solution  of  the  general  phase-space  ion 
distribution  were  proposed. 

The  above  work  was  considered  to  be  an  important  advance  in  the  field  by  a 
distinguished  expert  in  a  recent  review  of  the  subject  (cf.  D.  R.  Bates, 
Advances  in  Atomic  and  Molecular  Phvsics  (New  York,  Academic  Press,  1985),  vol . 


Theory  of  lon-Holecule  Collisions  at 
(1  eV-5  keV)/AMU 


by 

M.  R,  Flannery 


Published  as  a  chapter  in  the  book  "Swarms  of  Ions  and  Electrons  in 
Gases",  eds.  W.  Lindinger,  T.  0.  Mark  and  F.  Howorka  (Springer-Verlag,  1984) 
pages  103-125. 


Theory  of  Ion-Molecule  Collisions 
at  (1  eV-5  keV)/AMU 

M.  R.  Flannery 


School  of  Physics,  Georgia  Institute  of  Technology, 
Atlanta,  GA  30332,  U.S.A. 


;1.  Various  Aspects 

In  1on-ino1ecu1e  collisions  from  thermal  energies*  to  a  few  eV,  up  to  a  few 
keV* various  processes  such  as  chemical  reactions  (nuclear  rearrangements ) , 
i charge  transfer  (electronic  rearrangement)  and  rotational,  vibrational  and 
lelectronic  transitions  occur,  rarely  in  isolation  except  over  limited 
lenergy  ranges.  Inelastic  scattering  changes  the  internal  state  by  rotat- 
jional,  vibrational  and  electronic  excitation  of  one  or  both  collision 
iparteners  while  preserving  the  identity  of  each  partner.  In  contrast  to 
{reactive  scattering  (chemical  reactions,  three-body  Ion-molecule  assocla- 
Itlon,  charge  transfer,  dissociation.  Ionization)  which  changes  the  chemical 
{identities.  There  are  few  theories  which  elucidate  the  coupling  betiveen 
{reactive  and  Inelastic  mechanisms, since  development  of  feasible  theoretical 
treatments  of  each  of  the  orocesses  in  Isolation  oresents.  In  itself,  a 
considerable  challenge. 

The  electrostatic  interaction  V(R,r,Y)  between  an  ion  A'*’  at  relative  separa- 
i  tlon  ]g  from  the  center  of  mass  of  the  molecule  BC  with  intemuclear  vector 
ij|[,  depends  on  r  which,  in  general,  causes  vibrational  excitation  via  colll- 
islon  Induced  forces  along  on  the  angle  y  between  and  which  induces 
(via  a  torque)  rotational  excitation,  and  on  R  which  mainly  controls  the 
general  variation  of  the  corresponding  Inelastic  differential  cross  sections 
with  scattering  angle.  We  shall  see,  however,  that  vibrational  excitation 
can  be  strongly  enhanced  when  coupled  to  charge- transfer  channels.  Within 
various  regions  of  the  kinetic  energy  E  of  relative  motion  (or  center-of- 
mass  energy),  chemical  reactions,  charge  transfer,  rotational  and  vibration¬ 
al  excitation  may  be  all  or  partially  strongly  coupled,  or  else  be  all  effec¬ 
tively  decoupled  over  a  small  ('v-lOeV)  energy  range.  While  strono  couplings 
may  exist  only  within  all  reactive  mechanisms  (chemical  reaction  and  charge 


transfer  for  E  <  1  eV)  or  only  within  all  inelastic  mechanisms  (rotational 
and  vibrational  excitation}  in  ranee  1  <  E  (eV)  <  10-50,  i.e.,  essentially 
in  mutually  exclusive  energy  ranoes,  one  inelastic  mechanism  can  be  strongly 
coupled,  for  a  wide  range  of  interesting  cases,  with  only  one  strong  re¬ 
active  scattering  mechanism  (as  in  vibrational  excitation  caused  directly 
by  charge  transfer  for  10  <  E  (eV)  <  few  keV). 

A  plausible  upper  limit  to  reaction  at  kinetic  energy  E  of  relative  motion 
is  the  Langevin  cross  section,  Oj_(A  )  =  17o  '^[E(eV)]”  '2,  for  close 
spiralling  encounters  within  the  orbiting  radius,  '^<(4-6)  S  at  thermal 
energies,  centered  about  the  center  of  mass  of  the  molecule  of  averaged 
polarizability  a(X^).  At  thermal  energies,  the  cross  section  Oj.  for  pro¬ 
duction  [1]  of  Ar  in  (Ar^-H2)  collisions  is  Oq(^^)  200.  Also  HD^  and 

Dj  are  each  formed  in  merging  beam  (O^-HD)  collisions  [2]  with  O0(^^)  *^100, 
10  and  0.1  at  kinetic  energies  E  of  relative  motion  of  0.002,  1  and  8  eV, 
respectively.  The  E-variation  of  becomes  quite  precipitous  for  E  >  1  eV, 

w 

relative  to  the  much  slower  E-variation  of  Oj^.  Except  at  very  low  energies 

'E  <  1  eV,  charge  transfer  cross  sections  become  much  larger  than  a^, 

isince  attraction  can  occur  from  overlap  forces,  in  addition  to  the  weaker 
-4 

polarization  R  attraction  implicit  in  0|^. 

For  (H*,  D2)  (0^,  HD),  recrtion  and  charge- transfer  can  occur  for  E  (eV) 

|>  0.04  assisted  by  direct  rotational  and  possible  vibrational  excitation. 

I  For  (H*-H,(v=0,  J=0))  collisions,  in  contrast,  rotational  excitation  (with 

i 

Ithresholds  4E(J=0  ->■  2,4)  =  0.04,  0.14  eV)  and  vibrational  excitation  (with 

i 

jAE(v*0  -►  1,2,3)  =  0.516,  1.003,  1.461  eV)  are  alone  possible,  until  E 

jreaches  the  threshold  E^  1.83  eV  for  charge- transfer  and  reaction.  The 

icross  section  for  charge-transfer  (3,4]  then  slowly  rises  to  a  peak 

[■v  10  at  E  -v  1  keV.  Above  Eq  =  4.48  eV  direct  dissociation,  and  above 

|E^  11  eV  dissociation  via  electronic  excitation,  can  occur,  but  again 

•  '  + 

iboth  peak  in  the  keV  energy  region.  For  Li  -H2  collisions  at  E  <  E^  = 

110.2  eV,  only  rotational  and  vibrational  transitions  and  weak  direct 
idissociation  are  possible. 

Since  the  strong  valence  forces  in  Ht  tend  to  expand  the  H,  bond  length, 
vibrational  excitation  of  H2  by  H  is  much  stronger  than  in  Li  -H2  colli¬ 
sions  where  short  range  (forced  oscillator)  repulsive  forces  are  mainly 
responsible  for  the  small  vibrational  excitation  with  associated  excitation 
of  high  O'  levels.  Rotational  excitation  in  the  lighter  H^-H2  system  is, 
however,  weaker  by  comparison,  since  the  torque  is  less  and  the  triangular 
geometry  of  the  intermediate  complex  permits  less  flexibility  for  rota¬ 
tion  of  H,. 
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As  the  Ion  energy  or  kinetic  energy  of  relative  motion  E  increases  from 

%  several  eV,  rotational  excitation  decreases  since  the  rotational  period 
•11  -12 

*^(10  for  0-  -2  10  for  H-)s  becomes  much  greater  than  the  duration 
'V  R(A)  tH(A.M.U.)/2E(eV)l  '2  •'♦j;  of  collision  between  species 

of  reduced  mass  M(A.M.U.)  with  characteristic  length  R(%)  for  the  particular 
process.  The  cross  sections  ap  for  rotational  transitions  (in  Li^-Hg  colli¬ 
sions)  increase  from  threshold,  exhibit  a  peak<\<10  at  1-5  eV  and  then 
decrease  [5]  to-v-l  at  100  eV.  Since  the  vibrational  period  (10“^^ 

for  Brg  -  10"^^  for  Hgls,  for  vibrational  excitation  in  Li^-Hg  collisions, 
will  increase  from  thresholds  -vl  eV  mainly  via  v,J  -►  v',J'  transitions  say 
to  -vO.!  and  0.03  for  v  “  0  -►  1,2  respectively  [6]  at  E'vlOO-200  eV, 
and  then  decrease  thereafter.  The  H^-Hg  (v*0  -►  v')  system  involves  larger 
vibrational  cross  sections  with  maxima,  cr^  (%^)  6,  2.3,  and  0.8  for  v'  = 

<1,2  and  3,  respectively,  which  occur  [7]  at  E  ‘v  30  eV. 

Charge  transfer  collisions  [8]  between  similar  (AB^-AB)  systems  are  char¬ 
acterized  by  much  larger  cross  sections  which  decrease  only  very  slowly 

I  X  + 

as  E  is  raised  from  thermal  energies  to  a  few  keV;  for  N2-N2,  decreases 
![91  from  ^,70  at  E  -v-  1  eV  to  ''-25  at  a  few  keV.  Since  the  direct  and 
icharge  transfer  channels  can  be  in  near  energy-resonance  for  symmetric 
systems  and  are  closely  courle<l,  direct  vibrational  excitation  [8-10]  is 
[therefore  characterized  by  large  cross  sections,  ■v(10-20)  in  contrast 

ito  weaker  vibrational  excitation  via  forced  oscillator  collisions,  which  are 
characterized  by  relatively  smaller  cross  sections  -vio"^ 

The  key  difference  between  symmetric  (AB*-A8)  and  asymmetric  (A^-BC)  charge 

transfer  collisions  is  that  there  are  initially  two  (elastic)  paths  of 

approach  in  the  former  case,  and  charge  transfer  then  arises  from  the 

phase-change  caused  by  the  different  elastic  scattering  potentials  associa- 
I  ^ 

Ited  with  the  bonding  and  antibonding  states  of  the  [ABl,-complex  (states 

.which  reduce  to  the  gerade  and  ungerade  states  of  A2  in  A  -A  collisions). 

I Charge- transfer  in  the  asymmetric  (A^-BC)  case  involves  an  electronic 

'transition  between  states  with  different  asymptotic  energies  ('v-  1.83  eV 

for  H^-H2  and  11  eV  for  Li*-H2  collisions),  and  is  therefore  characterized 

by  a  threshold  dependent  cross  section  which  increases  to  somewhat  lower 

02  ^ 

cross  section  maxima  (o^  'v  10  A  at  E  1  keV  for  H  -H2  collisions  [3]). 

The  variation  with  E  of  many  ion-molecule  processes  appears  as  if  the  process 
were  energy  resonant.  This  is  because  the  manifold  of  vibrational  channels 
(in  charge- transfer,  for  example)  offers  a  band  of  pathways  which  are  in 
near-energy  resonance  either  with  the  incident  channel  or  with  a  particular 
product  (charge-transfer)  channel  which  is  strongly  coupled  to  the  incident 
channel. 


The  dominant  processes  for  unlike  systems  (A^-BC)  In  the  (1-100)  eV  range 
of  center-of-mass  energy  E  are  elastic  and  rotational  excitation,  while  for 
E<\.(1  eV  -  several  keV),  the  dominant  process  for  like  (AB^-AB),  apart  from 
'direct  elastic  scattering.  Is  charge- transfer  with  Its  associated  direct 
vibrational  excitation.  For  E  >  f ew  keV,  direct  electronic  excitation  to 
bound  electronic  and  to  dissociative  electronic  states  occurs  with  relative¬ 
ly  small  cross  sections  1  For  lower  E  <  1  eV,  chemical  reactions 
will  compete  quite  strongly  with  any  of  the  above  processes  when  operative. 

In  the  following  sections,  feasible  and  tractable  semiclassical  (classical 
path,  quantal  target)  theories  for  Inelastic  (rotational  excitation)  pro¬ 
cesses  alone,  and  for  reactive  processes  (synnetrlc  and  asymmetric  charge 
transfer)  coupled  to  direct  Inelastic  processes  (vibrational  excitation) 
will  be  outlined  and  discussed.  The  theories  are  particularly  appropriate 
for  the  energy  range  (1  eV  -  several  keV),  since  chemical  reactions  are 
specifically  excluded.  Full  quantal  treatments  are  not  presented  since 
their  application  Is  feasible  to  at  most  10  eV,  and  since  they  are  well 
^documented  elsewhere  [11-15].  Space  allows  only  the  mention  of  quasi- 
classlcal  trajectory  methods  [11,15,17],  valuable  for  highly  averaged 
quantities  In  rotational  excitation  In  non-reactIve  systems,  and  of  seml- 
;c1ass1ca1  S-matrIx  theory  [13,16]  and  classical  perturbatiur  theory  [17,18], 
which  are  all  based  on  classical  relative  and  Internal  motions  and  which 
are  valuable  towards  realization  of  the  conceptual  power  of  classical  mech¬ 
anics.  Apart  from  vibrational  excitation  In  colllnear  collisions,  the 
S-matrIx  theory,  however,  becomes  computationally  unfeasible  [11,18]  for 
three  dimensional  systems,  even  in  comparison  with  full  quantal  methods. 

2.  Rotational  and  Vibrational  Transitions  in  Non-Reactive  Collisions 
2.1  Quantal  Methods 

The  traditional  full  quantal  close  coupling  (CC)  treatment  of  rovibratlonal 
transitions  Is  the  generalization  [19]  of  the  formulation  [20]  of  Arthurs 
land  Dalgamo  for  scattering  of  a  rigid  rotor  (RR)  by  a  structureless  pro¬ 
jectile  to  vibrating  rotor  (VR)  collisions,  and  is,  by  now,  well  established. 
It  Is,  however,  practical  only  for  E  <  f ew  eV,  since  the  Inherent  partial 
wave  expansion  for  the  relative  motion  converges  only  very  slowly  as  E  Is 
Increased,  and  excessive  amounts  of  computer  time  are  required.  It  has,  as 
yet,  only  been  applied  [19]  to  rotational  and  vibrational  excitation  in  L1^ 
-H2  collisions  at  0.6  eV  and  1.2  eV  where  coupling  with  the  two  excited 
vibrational  channels  (v  *  1,2)  Is  extremely  weak,  and  vibrational  excitation 
Is  then  at  most  incidental. 


Many  simplified  quantal  schemes  (J-,-coupled  states  (CS)  [21],  infinite- 


order  sudden  (lOS)  [22],  effective  potentials  [23],  L-domlnant  [24])  based 
on  CC  have  therefore  been  Introduced  for  E  <  10  eV.  These  methods  seek  to 

*V 

reduce  the  (rotational)  dimensionality  of  CC  by  assuming,  for  example,  that 
the  Z-axIs  of  a  rotating  frame  always  lies  along  ^  (as  In  CS  [21]  and  In 
earlier  rotat1ng>atom  approximation  of  Bates  [25],where  transitions  are 
therefore  restricted  to  zero  change  aMj  In  the  Z-components  of  the  rota¬ 
tional  quantum  number  J)  or  by  assuming  that  the  orientation  ^  of  the  mole¬ 
cule  can  be  considered  fixed  during  the  collision  (as  In  the  sudden  or 
adiabatic-nuclei  approximation  [26]).  In  IDS,  a  combination  of  both  sudden 
and  CS  Is  adopted.  At  low  energies  (  <  few  eV),  vibration  may  either  be 
Ignored,  as  In  rigid  rotor  (RR)  approximations,  or  be  treated  exactly  (VR), 
within  the  above  methods  by  full  solution  of  the  vibratlonally  coupled 
equations. 

Since  these  dimensionality- reducing  (OR)  schemes  remain  based  on  a  partial- 
wave  analysis  they  are  also  very  time-consuming.  For  the  transition  J  -  0 
-*■  2  In  L1^-N2  collisions  at  4.23  eV,  which  Involves  close-coupling  of  many 
J-levels,  10^  partial  waves  were  required  In  a  recent  lOS  treatment  [27], 

I 

Under  the  added  requirement  that  reliable  Interaction  potentials.  In  general, 
ab-lnitio,  are  needed  as  Input,  It  Is  therefore  not  too  surprising  to  note 
ithat  only  a  few  special  systems  have  been  studied  In  detail  by  Uie  above 
jquantal  methods,  at  only  a  few  energies  chosen  to  coincide  with  measurement. 
For  example,  pure  rotational  and  rovibratlonal  transitions  have  been  studied 
In  the  collisions;  H^-Hg  by  CS-RR  [28]  at  3.7  eV,  by  lOS-VR  [29,30]  at 

4.67,  6  and  10  eV;  Ll'^-Hg  by  CC-RR  [31],  CC-VR  [19],  CS  [32]  for  E  5  1.2  eV 

and  by  lOS-VR  [33]  at  5.54  and  8.8  eV;  Ll'^-Ng  by  lOS-VR  [27]  at  4.23  eV. 

I 

Comparison  of  the  lOS-VR  cross  sections  [29]  with  rotational 1y  resolved 
measurements  [34]  for  H*-H2  (v*0,  J=0  -►  v* (0,1, 2, 3),  J')  Indicates  good  to 
{excellent  agreement  at  10  eV  for  the  variation  with  scattering  angle  of  the 
{vibrational  transition  probabilities  (ratio  of  cross  sections)  when  sutimed 
lover  final  levels  J',  and  excellent  agreement  for  the  relative  differential 
cross  sections  for  vibrational  excitation.  Only  scattered  agreement  [30] 

Is  noted,  however,  for  rotational  excitation  (0,0  0,J')  at  4.67  and  6  eV. 

The  agreement  between  lOS-VR  [33]  and  measured  [35]  vibrational  transition 
probabilities  (summed  over  O')  for  Ll^-Hg  collisions  Is  reasonable  only  at 
8.8  eV  and  very  poor  at  5.54  eV.  When  compared  with  the  only  measured 
absolute  Integral  cross  sections  (*^4.2  [5])  for  rotational  excitation 

In  L1*-H2  (0,0  -►  0,2)  collisions  at  12  eV,  (and  beyond),  the  lOS-VR  cross 
section  19.8  [33]  at  8.8  eV  Is  a  factor  of  5  higher;  and  the  lOS-VR 

cross  section  'v  0.3  [33]  for  vibrational  excitation  (v“0  -►  1)  at  8.8  eV 

Is  much  In  excess  of  the  measured  [6]  peak  values  ‘vO.08  which  are  spread 
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over  a  wide  range  around  150  eV. 

Many  of  the  above  reduced- dimensionality  quantal  schemes  have  varying 
degrees  of  success  e.g.,  the  most  accurate  of  these,  the  J^-coupled  states 
(CS)  method  [321  when  compared  with  CC  (191  predicts  for  Ll^-Hg  (J®0  -►  2) 
collisions  to  within  only  33%  at  0.6  -  1.2  eV,  (In  contrast  to  Its  success 
with  atom-molecule  collisions),  and  differential  cross  sections  In  sub¬ 
stantial  disagreement  at  Intermediate  angles.  Also  the  great  sensitivity 
[29,281  exhibited  by  any  method  to  the  choice  of  various  (reasonably  accu¬ 
rate)  Interactions  V(R,r,Y)  places  a  great  premium  for  highly  accurate  (time 
(Consuming)  calculations  of  these  anisotropic  surfaces  [19,27,291  over  an 
(extensive  range  of  (R,r,Y)  parameters. 

2.2  Semiclassical  Theory  (Elkonal  path-Quantal  target) 

There  is  therefore  some  need  for  more  accurate  yet  feasible  treatments  of 
the  energy  region,  1  <  E  (eV)  <  10,  and  great  need  for  treatments  feasible 
in  the  region  10  ev  <  E  <1  keV.  For  low  energies  <10  eV  only  relative 
{differential  cross  sections  have  been  measured  [34-361.  Some  measurements 
j(5,61  of  Integral  cross  sections  exist  for  rotational  and  vibrational 
excitation  In  Ll'*^-H2  collisions  at  center-of-mass  energies  E  In  the  range 
i(12-90  eV)  for  which  accurate  theory  has  not  yet  been  applied. 

'For  the  above  reasons,  and  In  order  to  explore  the  strengths  of  semlclass- 
!lca1  methods  which  naturally  focus  attention  on  colllslonal  time-scales 
land  which  naturally  provide  greater  physical  insight  to  various  levels  of 
approximation,  McCann  and  Flannery  [37-39]  have  developed  a  multistate 
orbital  (or  energy-conserving)  treatment  of  rotational  and  vibrational 
transitions  in  Ion-molecule  and  neutral -molecule  collisions  at  Intermediate 
and  high  Impact  energies.  Rather  than  seeking  schemes  which  reduce  the 
dimensionality  of  the  quantal  CC  within  a  phase  shift  analysis  of  the 
jrelatlve  motion,  the  basis  of  their  semiclassical  treatment  Is  Initially  a 
semiclassical  description  (via  an  elkonal  or  JWKB  wavefunctlon)  of  the 
relative  motion,  which  Is  then  coupled  to  the  full  quantal  response  of  the 
(Closely  coupled  Internal  modes  of  the  collision  partners,  via  the  Introduc¬ 
tion  at  each  R  of  an  Interaction  averaged  over  the  quantal  response  of  all 
of  the  Internal  modes.  Even  at  low  energies  ('vO.6  eV)  the  theory  realized 
considerable  success  [39]  when  compared  with  full  quantal  CC  differential 
and  Integral  cross  sections  [19], (Figure  1),  and  with  measured  relative 
differential  cross  sections.  The  procedure  apparently  Introduces  less  error 
than  does  the  best  (CS)  [32]  dimensionality-reducing  (DR)  quantal  scheme, 
as  evidenced  by  comparison  of  CS  [32]  and  EC  [39]  with  CC  [19]  for  L1^-H2 
at  0.6  eV.  Within  the  EC  method,  DR  approximations  may  In  addition  be 
made.  Figure  1  Indicates  the  excellent  agreement  between  the  semiclassical 
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Fig.  1.  Quanta!  [19,28]  and  semiclassical  [39]  differential  cross  sections 
ifor  elastic  (0  -►  0)  and  inelastic  (J=0  -►  2)  scattering  in  Li^-Hg  and 
collisions. 

1(39]  and  quantal  methods  [28]  under  the  common  CS  approximations  for 
collisions. 

The  method  is  extremely  efficient  and  accurate  in  that  transition  amplitudes 
iconverge  much  faster  in  the  impact-parameter  representation  than  they  do  in 
the  angular  momentum  (partial-wave)  representation.  It  is  therefore  quite 
a  natural  method  for  (1  eV  -  1  keV)  inelastic  collisions.  It  has  recently 
been  shown  [40,41]  to  be  also  highly  reliable  for  reactive  collisions  for 
E  <  10  eV. 

2.3  Hultistate  Orbital  (Energy  Conserving)  Method: 


In  contrast  to  customary  semiclassical  treatments  based  on  a  straight  line, 
V  the  classical  path  ^(t)  for  relative  motion  between  species  A  (ion  or  atom) 

and  BC  is  obtained  [37]  by  solving  Hamilton's  equations 

5  3Q^/at  •  3<H>/3P^  ;  3P^/3t  «  -  3<H>/3Q^  (1) 


directly  for  the  variation  with  time  t  of  the  generalized  coordinates 
Qj  (=  X,Y,Z)  for  Jg(t),the  vector  position  of  A  relative  to  the  center-of- 
mass  of  BC,  and  of  the  associated  generalized  momentum  PjCt)  for  motion 
under  the  (Ehrenfest)  averaged  Hamiltonian  [37], 

<H(Jg(t))>  -  Pj(J8(t))/2u  +  +  V(^,Jg(t))l'?(^.Jg(t))>,(2) 

where  u  Is  the  reduced  mass  of  the  (A«BC)  system.  The  system  wavefunctlon 
is  therefore  expanded  In  terms  of  the  complete  set  of  eigenfunctions 
of  the  unperturbed  Hamiltonian  H^  for  the  Internal  states  of  the  Isolated 
species  A  and  BC  of  energy  E^,  ®  ^n*n^*  ** 

*<<:•«)  ■  S  *„(«)  ♦nit)  e’‘PHS„<«)  -  »„(«>)  (5) 

n 

where  j|[  denotes  the  collective  Internal  coordinates  of  A  and  BC  In  the 
absence  of  Interaction,  and  where  separates,  via  the  Bom-Oppenhelmer 

approximation.  Into  associated  electronic,  vibrational  and  rotational  parts. 
The  relative  motion  Is  described  by  ^e  above  combination  of  functions  A^, 
and  which  possess  the  following  physical  characteristics.  The 
elkonal.  or  classical  action.  Is  the  solution  of  the  Hamilton- Jacobi  equa¬ 
tion  [37-39], 

'  '‘nl«>  ■  "S  •  »„«(«)  .  (*) 

dilch,  of  course.  Is  the  heavy-particle  limit  to  Schrodinger's  equation, 

[  ♦  (Zli/ti^)  V„„(J))*(!5)  >  ((VS„)^  -  1(v\))»(^)  5  (5) 

The  wavefunctlon,  i|;(R)  =  exp(iS„),  describes  relative  motion  In  the  static 
i channel  n  under  Interaction  V„„,  whidt  Is  a  diagonal  element  of 

Vij(ig)  -  <*i(^)|V{^.Jg)Uj{^)>  (6) 

iwhere  V(j(^,Jg)  is  the  Instantaneous  electrostatic  Interaction  between  A  and 

BC  at  separation  R.  The  wavenumber  and  kinetic  energy  of  relative  motion 

at  asymptotic  R  is  k„  and  1i‘^k„/2u  and  the  local  momentum  at  R  is  liK  (R). 

n  n  'x.  '\.n  'x* 

The  phase  X_  is  the  solution  of 
n 

^  •  W  ■  0 

which  represents  flux  conservation,  *  0),  In  static  channel  n,  since 

the  current  associated  with  the  exponential  factor  in  (3),  for  real 
and  X^,  Is 

ll,  •  5  (*/„)  .xp(-2X„)  ;js„  (8) 

when  (4)  Is  used  for  the  local  velocity  Thus,  flux  associated  with 
channel  n  of  (3)  Is  changed  only  via  transitions  to  other  states  f,  and  A^ 
Is  then  the  probability  amplitude  of  remaining  in  state  n.  On  assuming 
that  the  eikonal  S^(R)  contains  the  major  R-variation  of  (3),  the  stationary 
state  Schrodinger  equation  with  (3)  Inserted  then  yields  the  three- 


dimensional  equation  [37-39], 

5.  V*>  W«>  “p '<W  "'■'W  <’> 

nff 

The  first  term  in  (2)  is  the  kinetic  energy  of  relative  motion  and,  with  (3) 
Inserted,  the  second  term  reduces  under  the  assumption  that  the  orbits  in 
the  different  channels  do  not  differ  appreciably  (i.e.,  %  X^),  to  the 

r-averaged  interaction  [37-39], 

<*(«)>  ■  S  *  I  »J(!8)  »„<«)  Vnl*)  «p  (10) 

n  f 

The  time  t  (=  s/v(]g))  is  introduced  as  a  (durnny)  variable  to  measure  dist¬ 
ance  s  along  the  common  trajectory  i.e.,  »  K^aA^/as  =  (w/fi)  aA^/at, 

and  the  mean  momentum  is 

*  K^)  =  P(JJ)  =  uldR/dtl  (11) 

so  that  the  evolution  of  the  eikonal  difference  in  (9)  and  (10)  can  be 
written  as 

"  V.^(Jg(t))}  dt]/fi  (12) 

where  The  amplitudes 

C^(t)  »  A^(t)  exp[-(i/fi)  /  V^^(Jg(t))dt]  (13) 

to 

/<hen  ViSerted  in  (9),  therefore  satisfy  the  usual  set  of  coupled  equations, 

Ifi  ayat  •  S  C^(t)  V^^(Jg(t))  expde^^t/ti)  (14) 

n 

I  derived  for  the  internal  response  of  a  system  to  a  time  dependent  perturba- 
ition  V(;[|,Jg(t)) .  The  implicit  time  dependence  for  motion  under  the  averaged 
Hnteraction  (10),  which  reduces  with  the  aid  of  (12)  and  (13)  to, [37-39], 

'  9  ♦ 

!  <V(R(t))>  *  S  [|C^(t)r  +  S  C^(t)C^(t)V^^(JR(t))  exp(iE^^t/ti)  ]  (15) 

i  ^ 

'Is  determined  from  Hamilton's  equations  (1),  which  now  reduce  to  [37-39] 

;  aQ^/at  »  P^(t)/u  ;  aP^/at  =  -  S  s  C*(t)  C^(t)[aV^j^/aQ^]exp(ic^^t/ft) 

I  "  ^  (16) 

The  external  relative  motion  (16)  is  therefore  solved  self  consistently 
with  (14),  which  characterizes  the  quantal  internal  response  of  the  mole¬ 
cule,  at  all  times.  Two  important  consequences  [37,39]  of  this  procedure 
are  that  total  energy  is  conserved  at  all  times  durino  the  collision, 
(d<H>/dt=a<H>/at  *  0)  -  energy  is  therefore  being  continually  redistributed 
between  the  internal  modes  and  the  relative  motion  (a  valuable  asset)  -  and 
that  total  probability  S^|C^(t)|  is  naturally  conserved  at  each  point  on 
the  common  trajectory  provided  by  (15)  and  (16). 


The  differential  cross  section  per  steradian  is 


(17) 


■  ( k^k^ )  I  (2u/4i*2)  r)exD(  i  jc^-Jg)  |  V(^.Jg)  |  '('(jr.Jg)  ^ 

for  scattering  of  an  Incident  beam  directed  along  j(^(0,0)  Into  unit  solid 
angle  about  jc^(8.^).  When  the  main  contribution  to  (17)  arises  from  the 
stationary  phase  which  produces  the  classical  orbit  R(e,(»l,  then  the 
scattering  amplitude  and  hence  differential  cross  section  o^f(e,(^) 

Is  easily  deduced  from  flux  considerations,  as  [37,39] 

f^f(e.*)  =  11m  [A^(Jg)  exp  1S;()g)ItaJ,(e,«)]^^2  (18) 

ft  8^®® 

where  o”^(e,(>),  the  differential  cross  section  for  elastic  scattering  by 
^ff^V  channel  f,  or  by  <V{R)>  of  (15)  for  common  trajectories  JR^  ^(0»i>), 
can  be  determined  either  from  classical  or  quantal  procedures,  and  where 

s;<«)  •  ■  «>.•«  <”> 

is  the  classical  action  along  the  path  measured  relative  to  the  action  for 
undeflected  particles  with  the  same  wavenumber.  When  more  than  one  tra¬ 
jectory  contributes  to  a  given  scattering  angle  e,  then  the  primitive  result 
[37]  can  be  suitably  generalized  [39]  to  yield  (uniform)  differential  cross 
sections  which  exhibit  the  customary  oscillatory  structure  (Fig.  1)  due  to 

Interference  between  the  classical  actions  S'  (R)  associated  with  each 

n  'V* 

classical  path  which  contributes  to  a  specific  scattering  angle  e. 
Representative  comparison  with  full  CC  quantum  treatment  at  <1  eV  Is  shown 
In  Figure  1.  Implementation  to  higher  energies  (1  eV  -►  1  keV),  for  which 
the  above  treatment  was  mainly  designed  (since  full  CC  quantal  treatments 
then  become  prohibitive).  Is  underway. 

In  the  limit  of  high  Impact  energy,  the  scattering  becomes  Increasingly 
'concentrated  in  the  forward  direction,  and  the  classical  cross  section 
Og^(0.f)  ^f  used,  diverges  (albeit  extremely  slowly).  The  cross  section 
for  transitions  between  sublevels  and  associated  with 

'rotational  levels  and  J^,  respectively,  can  then  be  obtained  directly 
from  (17)  to  yield  [42,43], 

o^f(e,^)  »  (k^/k.)|k^  /\(K'p)[C^^'-^p,-)  -  «^^]pdp|2  (20) 

0 

In  the  heavy-parti cle/high-enerqy  limit,  where  is  the  spherical  Bessel 
function  of  order  a  =  (M^  -  M^),  where  K'^  =  (k^  -  k^)^  -  e^^/ti^v?,  and 
where  C^*‘(p,t)  are  the  solutions  of  (14)  evaluated  using  a  straight  line 
orbit,  and  trajectory, 

jj(p,^)  »  constant  ;  ^(t)  *  t  (21a) 

in  terms  of  the  Impact  parameter  |^.  The  Integral  cross  section  then  reduces 
to  [42,43] 
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g 


O^^(v^)  -  2Ti(yic^)  /  IC^^'-^Cp.-)  -  pdp  (21b) 

0 

Mhich  permits.  In  contrast  to  the  usual  expression  without  the  evalua¬ 
tion  also  of  elastic  scattering  cross  sections. 

When  the  coupling  between  the  classical  path  and  quantal  target  can  be 
Ignored  I.e. ,  when  energy  conservation  which  results  from  (15)  Is  therefore 

violated,  the  classical  trajectory  R(t),  e(t)  and  orbit  jR(e,^)  at  energy 

12  ^ 

E  «  ^  can  be  assumed  planar,  and  Is  determined  by  Integration  of, 

^M(dR/dt)^  »  E  -  (Vg(R)  +  EpWl  ;  de/dt  -  -pvJ/R^  (22) 

where  7^  Is  either  some  averaged  spherical  Interaction  Independent  of  the 
transition  amplitudes  C^(t),  or  Is  the  sfrfjerical  part  of  the  surface  V{r,R), 
so  that  #  remains  constant. 

The  above  energy  conserving  (EC)  classical  path  (CP)  method  reproduces  the 
rainbow  oscillatory  structure  and  Is  capable  of  very  effective  description 
of  Inelastic  rotational  rainbows  115]  characterized  as  structures  In  the 
final  rotational  distribution  of  differential  cross  sections  at  fixed 
angles  for  those  molecular  systems  (N^,  CO,  Na^)  associated  with  small 
rotational  constants  so  that  many  rotational  levels  must  be  Included. 

These  rotational  rainbows  arise  from  coherent  superposition  of  scattering 
from  different  orientation  angles  y.  As  the  final  level  Is  Increased 
the  rotational  rainbows  In  -►  transitions  shift  to  larger  angles  due 
to  hard  collisions  and  have  wider  widths  than  the  usual  diffraction  osc111a> 
tlons  which  occur  from  near  forward  scattering  by  the  repulsive  part  of  the 
'Interaction.  The  analysis  of  these  rainbows  by  the  pure  CC  quantal  method 
ils  not  at  all  feasible,  and  reduced  dimensionality  simplifications  as  CS 
land  lOS  may  be  too  Inaccurate. 

The  above  EC-CP  method  can  also  be  extended  [44]  to  cover  vibrational 
itransitlons  and  vibrational  rainbows.  For  vibrational  excitation  simpler 
time- dependent  semiclassical  schemes  based  on  an  uncoupled  classical  path 
(21a),  (22)  or  Its  three-dimensional  equivalent,  have  been  adopted  [45-48]. 
Billing  [49]  has  treated  vibrational  excitation  by  assuming  that  rotation 
and  relative  motion  can  be  determined  classically  from  various  Hamiltonians 
which  Include  an  asserted  average  over  the  one-dimensional  quantal  vibra¬ 
tional  response,  but  which  are  not,  however,  consistent  with  energy  con¬ 
servation. 


3.  Charge-Transfer  and  Associated  Vibrational  Excitation 

In  contrast  to  charge  transfer  in  (B^-C)  collisions  150]  involving  atomic 
species,  theoretical  development  of  viable  descriptions  of  AB  -AB,  A  -BC, 
BC*-A,  AB*-CD  charge-transfer  collisions  are  relatively  scarce  [8-10,51-61]. 
This  is  mainly  due  to  great  difficulty  of  accurate  determination  of  the 
various  electronic  interaction  matrix  elements  which  are,  in  general,  aniso¬ 
tropic  and  which  involve  orbitals  based  on  three  and  four  centers.  The 
added  complexity  introduced  by  acknowledgement  of  vibrational  and  even 
rotational  transitions  (which  proceed  via  the  anisotropic  part  of  the 
Interaction)  is  however  tractable,  and  given  the  electronic  interactions, 
semiclassical  multistate  calculations,  although  time-consuming,  are  feas¬ 
ible.  Bates  and  Reid  [8]  have  realized  the  first  tractable  and  simplified 
multi-vibrational  level  treatment  of  charge- transfer  between  like  systems 
(AB^-AB)  which  do  not  involve  an  electronic  transition.  Hedrick  et  al.[60] 
have  provided  the  first  detailed  account  of  A^-BC  and  BC  -A  charge- transfer 
collisions  which  involve  an  electronic  transition  and  have  shown  (as  with 
Bates  and  Reid  [8])  that  direct  vibrational  excitation  is  greatly  enhanced 
by  coupling  with  the  charge- transfer  channel. 


In  the  collision, 

A8^vj,j|)  +  AB(v^,J^)  AB*(vp^)  +  AB(vpJ^)  ;  (D)  (23a) 

AB(v^,Jp  +  AB*(vJ,J^)  ;  (X)  (23b) 

i  at  low  and  intermediate  E  <  f ew  keV,  it  can  be  assumed  that  only  the  ground 

^  + 

'electronic  states  of  each  species  participate, so  that  the  [AB]2  complex 

has  only  two  possible  electronic  states.  These  represent  the  active  elec- 

I  tron  remaining  either  attached  to  the  slow,  neutral,  target  AB  -  the  direct 

; (0)  channel,  (23a)  -  or  attached  to  the  (fast)  projectile  ion  -  the  exchange 

i channel  (X),  (23fa).  The  kinetic  energy  E  of  relative  motion  of  two  species 

|each  of  mass  M  (A.M.U.)  and  moving  with  relative  speed  u  (a.u.)  is  'v  25  M  u^ 

' (keV)  such  that  for  E  sufficiently  high  >  10  eV  when  the  collision  duration 

the  rotational  period,  electron- transfer  may  be  considered 

to  occur  at  a  fixed  orientation  ^  of  each  internuclear  axis,  so  that  for 

spherical  interactions,  rotation  of  the  target  and  projectile  remains 

unaffected  i.e.,  0^  ■  and  »  J^.  Vibrational  transitions  v|*^  -*•  v^*^ 

are  however  Important  since  the  vibrational  period. 

3.1  Two-State  Treatment  for  (A^-A)  and  AB*-AB  Charge-Transfer 

Of  all  the  charge  transfer  processes  involving  molecules,  (23)  is  closest 

in  spirit  to  the  symmetric  resonance  process  [50], 
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A*  (fast)  +  A  (slow)  -►  A^  (fast)  +  A  (slow)  ;  (D)  (24a) 

-  A  (fast)  +  A"^  (slow)  ;  (X)  (24b) 

between  atomic  species,  insofar  as  (23b)  proceeds  via  the  time-varyino  field 
generated  by  the  energy  splitting  e^(R)  -  e'(R)  associated  with  the  bonding 
and  antibonding  eigenfunctions,  the  electronic  states  of  the  [AB  -AB] 
complex,  which  become  symmetric  or  antisymmetric  with  respect  to  the  mid¬ 
point  of  the  separation  Jg  between  the  center  of  masses  of  AB^  and  AB,  only 
when  intemuclear  vectors  Rp  and  Rj  of  AB*  and  AB,  respectively,  are  equal 
and  parallel. 

Since  the  bonding-antibonding  splitting  becomes  naturally  decoupled,  via 
the  Bom-Oppenheimer  approximation,  from  the  Franck-Condon  factors  which 
couple  the  vibrational  states  of  ion  and  neutral,  (23a)  also  permits  feasible 
detailed  study  of  characteristics  of  convergence  in  the  number  of  vibra¬ 
tional  states  used  in  the  wavefunction  expansion.  Also,  direct  vibrational 
excitation  (23a)  is  achieved  via  its  close  coupling  with  the  charge- transfer 
channel  (23b)  normally  characterized  by  large  cross  sections  'v.  10  at 
'v  100  eV. 


For  straight  line  relative  motion,  the  cross  section  for  (24b),  which  lays 
the  conceptual  groundwork  fcr  i23b)  is  (50,8,9) 

®el  *  f 
0 

where  the  phase  difference  between  elastic  scattering  in  each  of  the  qerade 
and  ungerade  channels  (which  are  acknowledged  by  expanding  the  system  wave- 
function  as  a  linear  combination  of  gerade  (-•-)  and  ungerade  (-)  eigenfunc¬ 
tions  of  Ap  is  150,9] 

n*i  -  (l/2ti)  r(c*(R(t))  -  c-(R(t))]dt  (26) 

When  either  an  LCAO  (linear  combination  of  atomic  orbitals)  is  used  for  the 

gerade  and  ungerade  states  of  At,  or  when  the  system  wavefunctions  is  ex- 

n  X 

panded  in  terms  of  atomic  states  based  on  each  center,  then  (25) 

follows,  with  (50,9) 

nji  ■  (1/fi)  /"v(R(t))dt  (27) 

where  the  exchange  interaction  is  given,  in  general,  by  [50,9] 

(l-|<^fl^J>|2]V(R)  •  (<^J|V°|^°>  -  <4l4j><4f|V°kJ>l,  (28a) 

■  [<4?|v’'i4>  -  <4jl4f><4lv’^l4j>l.  (28b) 

The  initial  and  final  wavefunctions  for  the  active  electron  attached  to  the 
target  and  projectile  ionic  cores  are  ♦?();;^)  and  respectively,  with 


^2  •  15  <Cl*  ^  *  f  for  energy  resonance.  The  electrostatic  interaction 
for  the  direct  (0)  channels  associated  with  the  fast  ion  and  slow  target  is 
V*^(;^ll,Jg)  and,  for  the  exchange  (X)  channels  associated  with  a  fast  neutral 
and  slow  ion,  is  the  bra-ket  <|>  denotes  integration  over 

the  coordinate  (or  of  ff’e  active  electron,  the  integrations  in  (28) 
Involve,  in  general,  one-electron  two  center  integrals.  For  symmetric 
resonance  processes,  the  diagonal  terms  in  (28a)  and  (28b)  are  equal. 


By  comparison  of  (26)  and  (27),  the  gerade-ungerade  splitting  for  electronic 
states  in  the  LCAO  approximation  may  be  approximated  as 

e^R)  -  c"(R)  =  2V(R)  (29) 

On  assuming  in  (23)  that  only  one  pair  (vt  =  v^,  v^  =  v^)  of  vibrational 
levels  participate,  then  on  separating  the  wavefunction  for  the  isolated 
neutral  and  ionic  molecular  systems  into  their  electronic,  vibrational  and 
rotational  parts  (via  the  Born-Oppenheimer  approximation),  the  charge 

transfer  cross  section  for  (23b)  is  [52] 

•  - 

a*  »  2ir  /  sin^  pdp  (30) 

o 

where  is  the  diagonal  element  of  the  matrix  jP  with  element 

’’tn  •  =  'fn  ''nf  »’) 

where  each  (Franck-Condon)  positive  or  negative  amplitude  F.  is  the  vibra- 

+  in 

tional  overlap  for  the  AB  (V|)  -►  AB(Vj^)  transition.  The  electronic  inter¬ 
action  V(Jg)  in  (27)  now  involves,  in  general,  orbitals  on  four  centers. 

The  above  expression  (30),  derived  originally  by  Gurnee  and  Magee  [52],  is 
valid  in  the  low  energy  limit  only  when  the  ground  vibrational  levels  of 
both  ion  and  neutral  alone  participate  in  an  energy- resonant  process.  It 
Is  the  natural  generalization  of  the  electronic  result  (25)  for  this  case. 
For  vibrationally  excited  ions  and  ground  neutrals,  initially,  accidental 
degeneracy  can  occur  when  a  band  of  various  pairs  of  (vpV^)  vibrational 
states  is  in  close  energy- resonance  with  the  original  pair  (vT,  v.  =  0)  so 
that  (30)  is  inadequate,  even  at  low  energies  E.  Moreover,  as  E  is  raised, 
an  increasing  number  of  vibrational  levels  participate  in  the  processes  (23). 

3.2  Multilevel  Treatment  for  (AB^-AB)  Charge-Transfer 
The  wavefunction  for  the  time  dependent  response  of  the  internal  modes 
(represented  collectively  by  r)  of  the  (AB^-AB)  system  under  the  influence 
of  their  mutual  electrostatic  interaction  V(jr,JR(t))  is  expanded  as  [8] 

f(JC.t)  *  I  I  A;(t)  exp(-iE“t/ti)  (32) 

a*D,X  n  " 

In  terms  of  ♦^()[)  which  form  a  complete  set  of  molecular  eigenfunctions 
(with  electronic,  vibrational  and  rotational  parts)  of  the  unperturbed 


Hamiltonian  for  the  Isolated  molecular  systems  with  total  Internal 
energy  at  Infinite  center-of-mass  separation  R.  The  index  a  denotes 
nrhether  the  labelled  quantities  are  associated  with  the  direct  (D)  channels 
(23a)  or  with  the  charge-transfer  (X)  channels  (23b).  Insert  (32)  into  the 
time-dependent  Schrodinger  equation,  and  assume  that  electronically  aver¬ 
aged  quantities  such  as  the  overlap 

S(]g)  *  •'♦el  I  *61 V 

'''ox  D  X 

between  the  electronic  parts  i|(g|  of  the  wavef uncti ons  for  the  D  and  X 

channels,  and  the  diagonal  distortion  (secular)  terms 

are  all  independent  of  each  of  the  internal  nuclear  separations  igp  and 
of  the  projectile  P  and  target  molecules  T  i.e.,  S  and  ^  are  spherical 
in  R  so  that  rotational  state  of  P  and  T  remains  unchanged.  By  working  to 
lowest  order  in  S,  Bates  and  Reid  [81  then  found  that  the  amplitudes, 

C;(t)  =  Aj(t)  exp(i/R)  /  V°^^(jg(t))dt  (35) 

satisfied  the  following  set  of  coupled  equations  [8] 

iti  3Cj/n  »  V(R)  I  P^^  C“(t)  exp(ie^„t/ti)  (36) 

n 

where 

»(«)  ■  -  S(«)  (37) 

is  assumed  spherical,  where  is  the  product  (31)  of  Franck-Condon 
amplitudes  and  where  is  the  difference  between  the  com¬ 

bined  energies  of  the  vibrational  levels  of  the  ion  and  neutral  (the 
small  difference  between  the  rotational  energies  in  the  D  and  X  channels 
being  neglected).  The  label  a  in  (36)  is  D  when  a  is  X  and  vice  versa, 
i.e.,  the  0  channel  cannot  proceed  without  participation  of  the  X  channels 
and  vice  versa.  For  straight  line  relative  motion,  (21a)  along  the  Z-axis, 
the  cross  section  for  charge- transfer  (o  =  X) 

o®^(E)  =  2ir  /  |C“(p,t-»-)  |2  pdp  (38) 

reduces  to  (30)  for  the  energy  resonant  case,  (v^  =  v^,  v^  *  v^). 

Once  the  solutions  C®  are  obtained  from  (36)  subject  to  C?(p,t-*— )  =  5^., 
y  n  T  ni 

C^(p,t->-«)  *  0,  differential  cross  sections  are  given  by  (20)  with  a  =  0. 
The  use  of  (21b),  rather  than  (38),  permits  the  additional  evaluation  of 
direct  elastic  integral  cross  sections.  Also  the  assumption  of  straight- 
line  relative  motion  may  be  removed  by  the  orbital  (or  energy-conserving 
method)  of  McCann  and  Flannery  [37]  so  that  extension  to  the  lower  (eV) 
energies  can  be  obtained  [56,40].  It  is  worth  noting  that  the  above 
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treatment  ignores  electron  translational  factors  exp(iu’j[),  which  arise  from 
the  momentum  change  introduced  as  the  active  electron  jumps  from  the  target 
to  the  projectile  ion  and  which  become  important  mainly  in  the  keV  region. 

At  low  impact  speeds  u,  when  only  these  levels  n  which  form  a  band  in  near 
energy  resonance  with  the  initial  channel  need  be  included,  (since,  other¬ 
wise  the  exponent  in  (36)  oscillates  rapidly  to  yield  negligible  contribu¬ 
tion),  and  at  high  u,  exp  (ic^p  Z/tiu)  in  (36)  is  effectively  unity,  the  set 
(36)  may  then  be  solved  analytically  by  diagonal izati on  techniques  [8]. 


M- 


At  high  energies  this  procedure  yields  [8), 

®1f  “  '’fi  '^el^^^  •  V  “  °  ^ 

X 

such  that  the  charge  transfer  cross  section  a-x  when  summed  over  all  final 

X  '' 

vibrational  levels  f  is  simply  as  in  the  atomic  case  (25)  -  (27). 

Thus,  near  energy-resonant  channels  (c-^  %  0)  with  reasonable  vibrational 

X 

overlaps  control  the  process,  at  low  E,  while  a.^  at  large  E  is  mainly 
determined  by  ;  dependence  at  intermediate  E  results  from  coupling 
between  energy  defects  and  P^^. . 

The  main  limitation  to  application  of  this  theory  [(36)  and  (39)]is  the 
inherent  complexity  of  calculation  of  the  exchange  interaction  (37)  which 
involves  orbitals  based  on  four  centers  and  anisotropic  V.  Bates  and  Reid 
[81  extended  the  method  of  Firsov  for  the  splitting  (29)  and  calculated 
the  0  and  X  cross  sections  for  H2(v)  -  H2(v  =  0)  collisions  in  the  (C.M.) 
energy  range  1  eV  <  E  <  25  keV.  Flannery  and  associates  [9,10,51,53-56] 
.investigated,  in  detail,  N2“N2,  CO^-CO,  NO^-NO  and  O2-O2  collisions  by 
adopting  a  prescription  due  to  Sato  for  (29);  and  H2-H2,  D2-D2,  1^2”^2 
collisions  by  using  the  Bates-Reid  interaction.  Up  to  100  vibrational 
'Channels  may  be  coupled.  The  similar  vibrational  frequencies  of  the 
neutral  and  ionic  species  (N2,  CO)  results  in  a  grouping  of  the  energy 
defects  for  different  vibrational  channels  into  different  "bands"  separated 
by  approximately  one  vibrational  quanta.  The  largest  vibrational  overlaps 
F^j  are  in  the  near  resonant  band.  Dissimilar  vibrational  frequencies  for 
the  ionic  and  neutral  species  (NO,  O2)  prevent  this  grouping,  and  large 
vibrational  overlaps  are  spread  over  a  moderate  number  of  product  channels. 
For  (H2,  D2,  T2),  the  ionic  and  neutral  vibrational  frequencies  are  dis¬ 
similar  but  the  energy  defects  do,  however,  exhibit  band  structure  and  the 
vibrational  overlaps  as  spread  over  a  large  number  of  product  channels. 
These  characteristics  (energy  defects  versus  vibrational  overlaps)  of  these 
various  systems  are  reflected  in  the  variation  of  the  calculated  cross 
sections  [9,10,53-56]  with  energy  E. 


Snmp  calculated  cross  sections  [541  are  illustrated  in  Figure  2  for  the 
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first  ten  D  and  X  channels  with  the  smallest  energy  defects  for  *  0) 

^  ^2  ^^1  *  collisions.  At  low  E.  the  channels  in  close  energy  resonance 
have  the  largest  cross  sections  as  expected  from  (30).  Above  100  eV-ion 
energy,  ^e  vibrational  overlaps  start  to  play  an  important  role  and  the 
cross  sections  are  influenced  considerably  by  energy  defect  and  overlap, 
until  at  high  energies  5  keV),  the  vibrational  overlaps,  in  accord  with 
(39),  control  the  relative  cross  sections  which  Increase  with  increasing 
overlap  (Figure  2). 

Quantative  comparison  of  theory  and  measurements  is  difficult  in  that  the 
vibrational  population  of  the  reactant  ions,  produced  in  many  experiments 
by  high-energy  electron  impact,  can  be  questionable  due  to  competition 
between  direct  ionization  and  autoionization.  Figure  3  illustrates  compari¬ 
son  of  theoretical  charge  transfer  cross  sections  for  reactant  ions  hI 

(Vf  «  0,  1)  with  various  measurements  for  ions  produced  mainly  in  the  lower 
■  ^ 

levels  (v^  >  0,  1)  by  16  eV  electron  impact.  At  low  energies  ('v  1  ev) 
rearrangement  channels  which  produce  and  which  are  not  included  in  the 
above  theory,  seriously  compete  with  charge-transfer.  Comparison  between 
theory  and  measurements  for  other  systems  may  be  found  in  the  review  by 
Moran  [51]. 

Within  the  Bates-Reid  method,  Smirnov  and  associates  [57]  determined  the 
matrix  element  V(^)  of  (37)  by  asymptotic  methods  (which  are  based  on  the 
lassumption  that  electron  transfer  occurs  at  distances  R  large  compared 
with  molecular  dimensions).  They  were  then  able  to  extract  both  the 
ispherical  and  anisotropic  parts  of  V(Jg)  from  Hartree-Fock  one  electron 
iwavefuncti ons  for  a  wide  range  of  systems. 

iSInce  the  main  limitation  to  application  of  the  Bates-Reid  method  is  a  lack 
;Of  reliable  data  on  the  exchange  interaction,  Borkman  and  Cobb  [69]  have 
irecently  computed  ab-initio  potential -energy  curves  for  the  complex 

using  SCF  (self  consistent  field)  and  Cl  (configuration  interaction)  tech¬ 
niques.  This  is  the  only  comprehensive  study  yet  available  for  AB^-AB 
system,  is  extremely  time  consuming,  and  indicates  that  the  interaction 
V(]^)  is  extremely  anisotropic.  Rotational  excitation  is  therefore  expected 
to  be  important  for  E  <  100  eV  or  so.  Lee  and  De  Pristo  [70]  have  recently 
found  agreement  between  their  simple  model  of  V(]g)  and  the  ab-initio  results. 

Oe  Pristo  [40]  adopted  the  orbital  (or  energy-conserving)  method  of  McCann 
and  Flannery  (37,56],  rather  than  the  straight-line  relative  motion,within 
the  Bates-Reid  treatment,  and  has  shown  good-excellent  agreement  with  the 
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Fig.  2.  State- to-state  (vj  »  0,  *  0)  (v^,  v^)  charge-transfer  cross 

sections:  (v^,  v^j  label  each  curve  which  at  high  energy  become 
arranged  In  order  of  the  Indicated  vibrational  overlaps 
(from  [54]). 
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Fig.  3.  Present  theoretical  (labelled  v^  ■  0,1)  and  various  experimental 
(SH  [62],  K  (631,  MR  [64],  LBG  [65],  HA  [66],  •  [67],  (■,  v^  -  ( 
□  *  ■  1)  [68])  cross  sections  for  (X^r^,  v^)  -  HgCX^Eg,  v^ 

■  0)  charge- transfer  collisions  (from  [51]  and  [54]). 
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only  full  quantal  treatment  performed  by  Becker  [59]  for  O2-O2  charge  trans¬ 
fer  collisions  at  8  and  36  eV  (C.M.).  Finally  Sears  and  De  Pristo  [58] 
have  recently  developed  a  simple  number  scaling  relationship  In  which  they 
assert  that  the  entire  set  of  unknown  state- to-state  cross  sections  can  be 
generated  from  one  cross  section  and  two  parameters  based  on  overlaps, 
energy  defects  and  collision  durations.  It  essentially  represents  an 
effort  to  extend  the  exact  parametric  dependence  (39)  at  high  energies  down 
to  lower  energies  by  acknowledging  Interplay  between  and  P^. . 

3.3  Multlstate  Treatment  of  (A^-BC)  and  (BC*-A)  Charge  Transfer 
For  unlike  (asymnetrlc)  systems,  Hedrick  et  [60]  expand  the  system 
wavefunctlon  as  (32)  and  find  that  tte  amplitudes  A^*^(t)  satisfy  the  set 
of  coupled  equations, 

I  aJf,„  .  v“,  I  expdEj'St/*)  *  A* 


(40a) 


I  I  A;f„^  «p(tE“t/li)  *  V»,  A? 


where  F-  Is  the  Franck-Condon  amplitude  <vt|v  >  for  BC*(vt)  BC(v  ) 

aft  ft  ft  *11  •  n 

transitions,  and  where  confined  electronic  and  vibra¬ 

tional  energy  defects  between  the  0  and  X  channels.  The  matrix  elements 


(40b) 


(41a) 


-  [<1.x|V^*x>  ^  V*2(R)]  «n)t 
•  I<Fx|v'’|«>  i  V*,(l»l  F,„ 


(Alb) 


!  -  i<n\*>  3  S2,(r)i  f,„  6^^  ;  ^ 

represent  various  r-averages  of  the  direct  Interaction  Vn  between  A  and 

^  XV 

BC,  and  of  the  exchange  interactions  V  between  A  and  BC  over  post  and 

prior  electronic  wavef unctions  (♦  for  BC  alone,  ♦  for  BC*,  x  Tor  the  one- 
electron  orbital  of  A),  and  over  vibrational  and  rotational  wavefunctlons 
for  BC  and  BC*.  The  Incident  speed  as  before  (§  3)  Is  sufficiently  large 
that  rotational  transitions  are  neglected  (and  V^j  are  assumed  to  be  Inde¬ 
pendent  of  orientation  of  the  Internal  Intemuclear  vectors  of  BC  and  BC* 
and  hence  of  the  orientation  of  jg,  the  relative  channel  vector  between  A* 
and  BC).  In  contrast  to  the  elastic  bonding  and  antibonding  channels  In 
the  AB*-AB  case  (§  3.2),  charge- transfer  occurs  here  via  an  electronic 


"  •/  •/  V 


•  eHe  r9  1 


transition  between  distinct  electronic  states  1  and  2  of  different  asymptotic 
energies  and  represented  by  electronic  wavefunctlons  a  and  Also  f 
V22#  ^2^  f  vi|2»  contrast  with  the  synmetrlc  case  (§  3.2),  and  the  overlap 
Sio  cannot  In  general  be  neglected.  We  note,  with  the  aid  of  (34)  that  the 
set  (40)  reduces  to  (36)  for  equal  distortion  terms  Vjj*(=  ,  j  •  1.2) 

In  the  0  and  X  channels,  when  the  electronic  overlap  Sif  Is  neglected,  and 
when  Is  naturally  replaced  by  For  the  case  of  one  (Incident)  0 

channel  and  n-X  channels,  (40)  reduces  to  a  simplified  set  [60]. 

Hedrick  et  al.  [60]  calculated  differential  and  Integral  cross  sections 
(20)  and  (21)  from  the  solution  of  (40)  for  Ar^-Hg  collisions  at  Ion 
energies  between  0.9  and  3.5  keV.  The  Interactions  In  (41)  were  computed 
directly  from  simple  analytical  electronic  wavefunctions  for  H2  and  H2 
and  from  an  analytic  fit  to  the  numerical  3p-orbital  for  Ar.  Fair  agree¬ 
ment  with  experiment  was  obtained. 

More  recently,  Sidis  and  de  Bruijn  [61]  have  given  similar  theory  as  above 
but  used  asymptotic  methods  to  determine  the  various  Interactions  for  H2- 
Mg  collisions. 

4.  Conclusion 

He  have  therefore  discussed,  iri  s  1  of  this  Chapter  various  aspects  of  ion- 
molecule  collisions.  We  have  then  presented  in  5  2  a  theory  suitable  for 
differential  and  integral  cross  sections  for  rotational  transitions  with 
extension  to  vibrational  transitions.  In  !  3,  a  theory  which  couples 
reactive  (charge-transfer)  and  inelastic  (vibrational  excitation)  scattering 
!ls  described  for  charge- transfer  and  vibrational  excitation  in  like  (AB^-AB) 
iand  unlike  (A^-BC)  systems.  A  natural  contination  of  the  treatments  in  §  2 
'and  S  3  would  permit  the  additional  coupling  of  rotational  transitions  to 
the  present  theory  in  s  3.  This  coupling  would  be  mainly  Important  for 
collision  energies  E  <  100  eV,  but  would,  however,  place  great  emphasis  on 
iaccurate  determination  of  ion-molecule  bonding  and  antibonding  anisotropic 
Interactions,  which  at  this  stage  are  only  (partially)  known  for  the 
[29],  Li'^-Hg  [19],  Li'^-N2  [27]  and  H^-Hg  [69]  systems. 
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